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Abstract. A necessary and sufficient condition for a one-dimensional g-state 
n-input cellular automaton rule to be number-conserving is established. Two 
different forms of simpler and more visual representations of these rules are 
given, and their flow diagrams are determined. Various examples are presented 
and applications to car traffic are indicated. Two nontrivial three-state three- 
input self-conjugate rules have been found. They can be used to model the 
dynamics of random walkers. 



1. Introduction 

A one-dimensional cellular automaton (CA) is a discrete dynamical system, 
which may be defined as follows. Let s:ZxN^Qbea function that satisfies 
the equation 

s(i,t+l) = f(s(i-re,t),s(i-r e + l,t),... ,s(i + r r ,t)), (1) 

for all i £ Z, and all t £ N. Z is the set of all integers, N the set of nonnegative 
integers, and Q a finite set of states, usually equal to {0, 1, 2, . . . ,q— 1}. s(i, t) 
represents the state of site i at time t, and the mapping / : Q r *+'>+ 1 _ > Q is 
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the CA evolution rule. The positive integers re and r r are, respectively, the left 
and right radii of the rule. In what follows, / will be referred to as an n-input 
rule, where n is the number re + r r + 1 of arguments of /. Following Wolfram 
pf|; to each rule / we assign a rule number N(f ) such that 



(Sl,X2,... ,x n )£Q 



N(f)= V f(xi,x 2 ,... ,x n )q^~ 



Cellular automata (CAs) have been widely used to model complex systems in 
which the local character of the evolution rule plays an essential role [§, [|, f|, ||] . 
In this paper, we will only consider finite CAs, and replace the set Z by the set 
Zl of integers modulo L. Any element of the set Q L will be called an L- cycle 
or a cyclic configuration of length L. 

Recently, following Nagel and Schreckenberg Q , many authors proposed vari- 
ous CA models of highway traffic flow. In the simplest models of this type, 1l 
represents a one-lane circular highway, and Q is equal to {0, 1}. According to 
the value of s(i, t), it is said that, at time t, site i is either empty, if s(i, t) = 0, 
or occupied by a car, if s(i, t) = 1. In order to complete the description of the 
model, an evolution rule has to be defined. For the sake of simplicity, we might 
assume that, at each time step, all cars move to the right neighboring site if, 
and only if, this site is empty. It is not difficult to verify that this rule coincides 
with Rule 184 defined by 

/i 84 (0,0,0) = 0, / 184 (0,0,1) = 0, / 184 (0,1,0)=0, / 184 (0,1,1) = 1, 
/i 84 (l,0,0) = l, / 184 (1,0,1) = 1, / 184 (1,1,0) = 0, / 184 (1,1,1) = 1. 

This traffic rule does not allow cars to enter or exit the highway. Therefore, 
during the evolution, the number of cars should remain constant, that is, start- 
ing from any initial cyclic configuration of length L, and for all teM, Rule 184 
should satisfy the condition 

/i84 (s(l, t), s{2, t), s(3, t)) + /i 84 (s(2, t), s(3, t), s(4, *))+••• 

+ /i 84 (s(L, t), s(l, t), s(2, t)) = s(l, t) + s{2, *) + ••• + s{L, t). (2) 

Since Rule 184 may also be written 

fim{xi,X2,x 3 ) = x 2 + min{a;i, 1 - x 2 } - mm{x 2 , 1 - x 3 }, (3) 

condition (||) is clearly verified. If we had assumed that cars, instead of moving 
to the right, had to move to the left, we would have found that the evolution 
rule would have been Rule 226. The relations 

f 22 e(xi,x 2 ,x 3 ) = fis4(x 3 ,x 2 ,x 1 ) 
f 22 6(xi,x 2l x 3 ) = 1 - /i 84 (l - xi, l—x 2 ,l- x 3 ), 

show that these two rules are not fundamentally distinct. 

Rules 184 and 226 are the simplest nontrivial examples of number-conserving 
CAs. In a previous paper 0, in order to find all the deterministic car traffic 
rules, we determined, up to n = 5, all the 2-state number-conserving CA rules. 
This was done using a different technique than the one presented in this paper. 
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For n > 3, not all these rules did correspond to realistic car traffic rules since 
they allowed vehicles to move in both directions. We found then that it was 
more appropriate to interpret all these 2-state number-conserving CA rules as 
evolution operators of one-dimensional systems of distinguishable particles. 
The purpose of this paper is to derive a necessary and sufficient condition for 
a one-dimensional g-state n-input CA rule to be number-conserving. We will 
then give examples of such rules and indicate some of their applications. Our 
result is an illustration of a general theorem on additive conserved quantities 
established by Hattori and Takesue ||. 

2. Number-conserving rules 

Definition 2.1. A one-dimensional g-state n-input CA rule / is number-con- 
serving if, for all cyclic configurations of length L > n, it satisfies 

f(xi,x 2 ,... ,x n -i,x n ) + f(x2,x 3) ... ,x n ,x n+ i)-\ 

+ f(%L, xi ... , x„„ 2 , ar n -i) = xi + x 2 H h x L . (4) 



Theorem 2.1. A one- dimensional q-state n-input CA rule f is number-con- 
serving if, and only if, for all (x\, x 2 , . . . ,x n ) € Q n , it satisfies 

n-l 

f{xi,x 2 , ... ,X„) = Xi + ^ (/( 0,0, ■ ■ ■ ,0 ,X 2 ,X 3 , . . . ,X n - k +l) 

fe=1 ' k ' 

-/(0, 0, . . . , 0, X\, X 2 , . . . , Xn-k)) ■ (5) 

V ' 

k 

To simplify the proof we will need the following lemma. 

Lemma 2.1. If f is a number-conserving rule, then 

/(0,0,... ,0)=0. (6) 

Write Condition (^) for a cyclic configuration of length L > n whose all elements 
are equal to zero. □ 

To prove that Condition (|^) is necessary, consider a cyclic configuration of 
length L > 2n — 1 which is the concatenation of a sequence (xi,x 2 , ... , x n ) 
and a sequence of L — n zeros, and express that the n-input rule / is number- 
conserving. We obtain 

/(0,0,... ,0,ari) + /(0,0,... , 0, x u x 2 ) + ■ ■ ■ 

+ f(xi,x 2 , ... ,x n ) + f(x 2 ,x 3 , ... , x n , 0) H 

+ f(x n ,0, . . . ,0) = xi + x 2 H \-x n , (7) 
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where all the terms of the form /(0, 0, . . . ,0), which are equal to zero according 
to (jfy, have not been written. Replacing x\ by in (Q) gives 

/(0,0,... ,0,x 2 ) + --- + f(0,x 2 ,... ,x n ) 

+ f(x 2 ,x 3 , . . . ,x n ,Q) H h f(x n ,0, ... ,0) 

= x 2 H Yx n . (8) 

Subtracting (§) from §) yields ©. 

Condition (^|) is obviously sufficient since, when summed on a cyclic configu- 
ration, all the left-hand side terms except the first cancel. □ 

Remark 2.1. The above proof shows that if we can verify that a CA rule / 
is number-conserving for all cyclic configurations of length 2n — 1, then it is 
number-conserving for all cyclic configurations of length L > In — 1 . 

The following corollaries are simple necessary conditions for a CA rule to be 
number-conserving. 

Corollary 2.1. If f is a one- dimensional q-state n-input number- conserving 
CA rule, then, for all x G Q, 

f(x, x, . . . , x) = x. (9) 
To prove (||), which is a generalization of (||), write Condition ([|) for x\ — x 2 = 

X fi X. I — I 

Corollary 2.2. If f is a one- dimensional q-state n-input number- conserving 
CA rule, then, 

J2 f(xi,x 2 ,... ,x n ) = ±(q-l)q n . (10) 

(x 1 ,x 2 ,... ,i„)£2™ 



When we the sum over (x\, x 2 , ■ ■ ■ ,x n ) G Q n , all the left-hand side terms 
except the first cancel, and the sum over the remaining term is equal to (0 + 
l + 2 + --- + (q-l))q n - 1 = \{q-l)q n . ' □ 

It is possible to give an interesting alternative proof of Relation ( |To| ) . Consider 
a De Bruijn cycle of length q n . Such a cycle contains all the different n-tuples 
{x\,x 2 ,... ,x n ) G Q™. The existence of De Bruijn cycles is related to the 
existence of Eulerian circuits on a De Bruijn graph G n -i- Such a graph has 
vertices and q n arcs. The number of arcs leaving a vertex is the same as 
the number which arrive. Each vertex is labeled by an (n — l)-tuple over Q, and 
the arc joining the vertex (x\, x 2 , ■ ■ ■ , x n -i) to the vertex (x 2 , . . . , x n -i, x n ) is 
labeled (x\,x 2 , . . . ,x n -i,x n ). For example, 



(0,0,0,0,1,0,1,0,0,1,1,0,1,1,1,1) 
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is a De Bruijn cycle for q — 2 and n — 4. It can be shown that there exist 
q~ n (ql) q distinct De Bruijn cycles jj). Since all the elements of Q appear 
an equal number of times in a De Bruijn cycle, Condition (0) implies 

f(x!,X 2 ,... ,X n ) = (0 + l + 2 + --- + (q-l))q n -\ 

and (0) follows. □ 



3. Examples of number-conserving rules 

Using Condition (^J) we can determine all the number-conserving CA rules for 
fixed values of q and n. These rules will be considered as operators governing 
the discrete dynamics of systems of particles whose total number is conserved. 
The particles occupy the cells of a one-dimensional periodic lattice subject to 
the condition that, at a given time, no more than q — 1 particles may occupy 
the same cell. Among these rules some have similar properties, which can 
be exhibited using the operators of reflection and conjugation. These two 
operators denoted, respectively by R and C, are defined on the set of all one- 
dimcnsional g-state n-input cellular automaton rules by 

Rf{xi,x 2 , ■ ■ ■ ,x n ) = f(x n ,x n -\, ... ,Xi) 

C f(xi,x 2 , ... ,x n ) = q- 1 - f(q - 1 - Xi,q - 1 - x 2 ,.. . ,q-l- x n ). 

It is clear that, if / is number-conserving, then Rf, C f, and RC f = CRf 
have the same property. Rules / and R f govern identical dynamics, the only 
difference is that, if for one rule particles flow to the right, for the other rule, 
they will flow in the opposite direction. To understand the difference between 
rules / and C / , let us assume that a cell, which contains k particles (1 < k < 
q — 1) contains q—l — k holes. Then, conjugation may be viewed as exchanging 
the roles of particles and holes. That is, if / describes a specific motion of 
particles then C f describes the same rule, but for the motion of holes. 
When the number of states and number of inputs are not very small, the 
dynamics of the particles is not clearly exhibited by the rule table of a number- 
conserving CA rule. A simpler and more visual picture of the rule is given 
by its motion representation. Such a motion representation may be defined as 
follows. List all the neighborhoods of a given occupied site represented by its 
site value s. Then, for each neighborhood, indicate the displacements of the s 
particles by a nondecreasing sequence of s integers (v\, v 2 , ■ ■ ■ ,v s ) representing 
the different velocities of the s particles. Velocities are positive if particles move 
to the right and negative if they move to the left. Alternatively, to the sequence 
(vi, «2, ■ ■ ■ ,v s ), we can substitute arrow(s) joining the site where are initially 
located the s particles to the final positions of the particle(s). A number above 
the arrow indicates how many particles are moving to this final position. To 
simplify a bit more these representations, we only list neighborhoods for which, 
at least one velocity Vj (j = 1,2,... , s) is different from zero. For example, 
these two forms of the motion representation of 2-state 3-input Rule 184 are 



i 

10 (1) and 10. 
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Since, for Rule 184, particles move only to the right, there is no need to indicate 
the state of the left neighboring site of the particle. Many other examples of 
motion representations for 2-state 4- and 5-input rules are given in 

Remark 3.1. In many applications, as one-way road car traffic, which pro- 
hibits passing, we have to assume that the particles are distinguishable. We 
have, therefore, to label them using an increasing sequence of integers, in order 
to be able to follow the motion of each individual particles. Instead of describ- 
ing formally the labeling process, we will give a simple example. Suppose we 
have the configuration: 










1 





2 





2 







We first replace each particle in 


each 


occupied cell by the s 








• 




• • 




• • 






and then label consecutively all the • to obtain: 








1 




2 3 




4 5 







ymbol • 



If now, we assume that the motion representation is 

01 (-1) 020 (-1,1) 021 (-1,1) 120 (0,1) 121 (0,1) 

220 (0,1) 221 (0,1) 022 (-1,0) 022 (-1,0), 
then, applying this rule, we obtain the new labeled configuration: 





1 




2 




3 4 




5 





which corresponds to the new configuration: 





1 





1 





2 





1 





Note that the above labeling convention assumes that particles cannot jump 
above each other, and, therefore, the ordering of labels remains unchanged 
after each iteration. We use this convention to ensure uniqueness of the motion 
representation for number-conserving CA rules. 

3.1. Three-state two-input number-conserving rules 

There are only 4 three-state two-input number-conserving rules, which are 
listed in Table 1. 

For radii pair (rg,r r ), we have two possible choices, either (0, 1) or (1,0). In 
each case, we have chosen the pair simplifying most the motion representation. 
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rule number 


base 3 representation 


(rt,r r ) 


19305 


222111000 


(0,1) 


15897 


210210210 


(1,0) 


18561 


221110110 


(1,0) 


16641 


211211100 


(0,1) 



Table 1. Three-state two- input minimal number-conserving rules. 

Rules 19305 and 15897 coincide with the identity. The motion representations 
of Rules 18561 and 16641, which are obtained from one another by reflection, 
are, respectively, 

20 (0,1), 21 (0,1) and 02 (-1,0), 12 (-1,0), 

or 

ii ii 
20, 21, and 02, 12. 

For both rules, if a site is occupied by 2 particles, one of them move to a 
neighboring site, except if this site is already occupied by 2 particles. 

3.2. Three-state three-input number-conserving rules 

There are 144 three-state three-input number-conserving rules. They can be 
divided into 48 equivalence classes under the dihedral group generated by the 
operators R and C. Each class will be represented by the rule having the 
smallest rule number, and which will be called the minimal rule of the class. 
Tables 2 and 3 list rule numbers and corresponding base-3 representations of 
all 48 minimal rules. 

Among these minimal rules, 18 eventually emulate the identity. Their reference 
numbers are: 27, 28, 29, 30, 31, 32, 33, 35, 37, 39, 40, 42, 43, 44, 45, 46, 47, 
48. The motion representations of the remaining 30 rules are given in Table 4. 
In order to have a global view of the properties of the various minimal rules, 
Figures [l] and ^ represent the flow diagrams of the 30 minimal rules whose 
motion representations are listed in Table 4. If p is the average particles density 
and i> av the average particles velocity, a flow diagram shows how the flow /cw av 
varies as a function of p. They are common practice in car traffic theory. 
Concerning these flow diagrams, one point has to be stressed. For q > 2, that 
is, if more than one particle can occupy a site, it is not always possible to 
choose the left and right radii of a rule so as to have u av = for p = 1. This 
is due to the fact that, choosing (rg,r r ) = (1, 1), for some rules, pv av = —0.5 
when p = 1. Any other choice will either increase or decrease i> av by one unit. 
The flow of the rules which emulate the identity, is equal to zero for all values 
of p. 

Remark 3.2. All the rules for which the flow is non- negative for all values of 
the average density, could be a priori considered as deterministic one-lane high- 
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£ 1 

rcicrcncc number 


rule number 


base-3 representation 


1 


6159136430181 


210210210210210210210210210 


2 


6159523870341 


210210211210210211210210100 


3 


6169984499133 


210211211210211100210211100 


4 


6171146819301 


210211221210211110210100110 


5 


6171534259461 


210211222210211111210100000 


6 


6181994888253 


210212222210212000210101000 


7 


6201360244413 


210221211210110100210221100 


8 


6202522564581 


210221221210110110210110110 


9 


6202910004741 


210221222210110111210110000 


10 


6213370633533 


210222222210111000210111000 


11 


6436931290101 


211210100211210211211210100 


12 


6447391918893 


211211100211211100211211100 


13 


6448554239061 


211211110211211110211100110 


14 


6448941679221 


211211111211211111211100000 


15 


6450878870973 


211211200211211200100211200 


16 


6452428631301 


211211211211211211100100100 


17 


6459402308013 


211212111211212000211101000 


18 


6462889260093 


211212211211212100100101100 


19 


6478767664173 


211221100211110100211221100 


20 


6480317424501 


211221111211110111211110000 


21 


6482254616253 


211221200211110200100221200 


22 


6483804376581 


211221211211110211100110100 


23 


6490778053293 


211222111211111000211111000 


24 


6491940373461 


211222121211111010211000010 



Table 2. Three-state three-input minimal number-conserving rules, Part 1. 
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reference number 


rule number 


base-3 representation 


25 


6494265005373 


211222211211111100100111100 


26 


6495427325541 


211222221211111110100000110 


27 


6726349098981 


212211000212211111212100000 


28 


6729836051061 


212211100212211211101100100 


29 


6736809727773 


212212000212212000212101000 


30 


6740296679853 


212212100212212100101101100 


31 


6757724844261 


212221000212110111212110000 


32 


6761211796341 


212221100212110211101110100 


33 


6768185473053 


212222000212111000212111000 


34 


6769347793221 


212222010212111010212000010 


35 


6771672425133 


212222100212111100101111100 


36 


6881331565845 


220100211220211100220211100 


37 


6893341954965 


220101222220212000220101000 


38 


6912707311125 


220110211220110100220221100 


39 


6924717700245 


220111222220111000220111000 


40 


6956093445525 


220121222220010000220121000 


41 


7158738985605 


221100100221211100221211100 


42 


7170749374725 


221101111221212000221101000 


43 


7174236326805 


221101211221212100110101100 


44 


7202125120005 


221111111221111000221111000 


45 


7205612072085 


221111211221111100110111100 


46 


7233500865285 


221121111221010000221121000 


47 


7448156794485 


222101000222212000222101000 


48 


7479532539765 


222111000222111000222111000 



Table 3. Three-state three-input minimal number-conserving rules, Part 2. 
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reference number 


motion representation 


1 


1 (-1) 2 (-1,-1) 


2 


10 (1) 20 (-1,1) 21 (-1,0) 22 (-1,0) 


3 


20 (-1,1) 21 (-1,0) 22 (-1,0) 


4 


01 (-1) 20 (-1,1) 21 (-1,0) 22 (-1,0) 

V/ V'/ V'/ V'/ 


5 


10 (1) 11 (1) 20 (1,1) 21 (0,1) 


6 


20 (1,1) 11 (1) 21 (0,1) 


7 


11 (-1) 20 (-1,1) 21 (-1,0) 22 (-1,0) 


8 


01 (-1) 11 (-1) 20 (-1,1) 21 (-1,0) 22 (-1,0) 


9 


10 (1) 20 (1,1) 21 (0,1) 


10 


20 (1,1) 21 (0,1) 


11 


10 (1) 2 (-1,0) 


12 


2 (-1,0) 


13 


01 (-1) 2 (-1,0) 


14 


10 (1) 11 (1) 20 (0,1) 21 (0,1) 

V/ V/ V"/ V"/ 


15 


02 (-1,-1) 12 (-1,0) 22 (-1,0) 


16 


10 (1) 11 (1) 020 (-1,1) 021 (-1,1) 120 (0,1) 




121 (0,1) 220 (0,1) 221 (0,1) 022 (-1,0) 


17 


11 (1) 20 (0,1) 21 (0,1) 


18 


020 (-1,1) 021 (-1,1) 120 (0,1) 121 (0,1) 




220 (0,1) 221 (0,1) 11 (1) 022 (-1,0) 


19 


11 (-1) 2 (-1,0) 


20 


10 (1) 20 (0,1) 21 (0,1) 


21 


11 (-1) 02 (-1,-1) 12 (-1,0) 22 (-1,0) 


22 


10 (1) 020 (-1,1) 021 (-1,1) 120 (0,1) 121 (0,1) 




220 (0,1) 221 (0,1) 022 (-1,0) 


23 


20 (0,1) 21 (0,1) 


24 


01 (-1) 20 (0,1) 21 (0,1) 


25 


020 (-1,1) 021 (-1,1) 120 (0,1) 121 (0,1) 




220 (0,1) 221 (0,1) 022 (-1,0) 


26 


01 (-1) 020 (-1,1) 021 (-1,1) 120 (0,1) 121 (0,1) 




220 (0,1) 221 (0,1) 022 (-1,0) 122 (-1,0) 


34 


01 (-1) 21 (0,1) 


36 


21 (-1) 20 (-1,1) 21 (-1,0) 22 (-1,0) 


38 


11 (-1) 21 (-1) 20 (-1,1) 21 (-1,0) 22 (-1,0) 


41 


21 (-1) 2 (-1,0) 



Table 4. Motion representations of the 30 minimal which do not emulate the identity. 
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Figure 1. Flow diagrams, Part 1. 



way car traffic. Consider for instance Rule 6171534259461 (reference number 
5), the second form of its motion representation is 

112 1 

rv rv rv rv 

10 11 20 21. 

That is, as many particles as possible move from one site to the right neigh- 
boring site. This rule, which is the simplest three-state rule generalizing Rule 
184, may be written 

f(xi, x 2 ,x 3 ) = x 2 + min{a;i,2 - x 2 } - min{a;2, 2 - x 3 }. 

This expression suggests a further generalization. The g-state rule defined by, 
for all (x 1 ,x 2 ,x 3 ) G {0, 1, . . . ,q — l} 3 , 

f(x 1 ,x 2 ,x 3 ) = x 2 + mm{xi,q - 1 - x 2 } - min{x 2 ,q - 1 - x 3 }, 

can be viewed as the following car traffic rule: Each cell represents a section 
of a one-way road between two traffic lights. The number of states q measures 
the maximum capacity of that section. The above evolution rule describes the 
way cars move from one section to the next. As for elementary CA Rule 184, 
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Figure 2. Flow diagrams, Part 2. 



the flow is maximal for p = 0.5. Of course, and this goes beyond our purpose 
here, we could, as in the Nagel-Schreckenberg model ||, introduce some noise 
and say that some cars, which could move to the next section, do not do so 
with a probability p. The existence of traffic lights could be used to monitor 
traffic in order to increase the flow, and, with this in mind, we could extend to 
<7-state rules classes of models we studied in a recent paper [n0| . 



Remark 3.3. As mentioned above, conjugation exchanges the roles of parti- 
cles and holes. Therefore, if a rule / is self-conjugate, we can always choose 
the rule radii (ri, r r ) such that the point (0.5,0) is a center of symmetry of its 
flow diagram. Rules whose reference numbers are 26 and 34 are self-conjugate 
and their flow diagrams have this property. In Q we studied a few 2-state 
self-conjugate number-conserving rules. Some of these rules, which allow mo- 
tion in both directions, govern the dynamics of ensembles of one-dimensional 
pseudo-random walkers. Starting from a random initial configuration, whose 
density p is exactly equal to 0.5, we followed, in the limit set, which is reached 
after a maximum of L/2 iterations, the position of a specific particle as a func- 
tion of time on a ring of length L — 5000 for 5000 time steps. The walks are 
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represented in Figures Since L is finite, the random walks are periodic 

in time. The first walk has a period equal to L, while the period of the second 
one is 2L. Note that these walks are deterministic, their random character 
comes from the randomness of the initial configuration. 

















\ M 
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1 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
t 

Figure 3. Random walker evolving according to Rule 6495427325541 (no 26) 
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Figure 4. Random walker evolving according to Rule 6769347793221 (no 34) 



4. Conclusion 

We have established a necessary and sufficient condition for a q-state n-input 
CA rule to be number-conserving. For given values of q and n, this result 
allows to find all the rules possessing this property. As an example, we have 
determined all the three-state three-input number-conserving CA rules. We 
have listed their motion representation and studied their flow diagrams. All 
the rules for which the flow is non-negative for all values of the average density, 
can be considered as deterministic one-way road car traffic rules in which the 
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cells represent road sections whose car capacity is equal to q — 1. Among the 
148 three-state three-input rules there exist two nontrivial self-conjugate rules 
which mimic the dynamics of an ensemble of random walkers. 
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